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7 ABSTRACl^ 

control strategy for unattended wind turbine operation is a potentially 
important aspect of wind energy production that has thus far escaped treatment in the 
literature. Early experience In automatic operation of the Sandia 17-m VAWT has 
demonstrated the need for a systematic study of control algorithms. To this end, a 
computer model has been developed using actual wind time series and turbine performance 
data to simulate the power produced by the Sandia 17-m VAWT operating in automatic con- 
trol- TfiF mddeT Tia"s been used to investigate the influence of starting algorithms 
on annual energy production. The results indicate that, depending on turbine and 
local wind characteristics, a bad choice of a control algorithm can significantly re- 
duce overall energy production. The model can be used to select control algorithms 
and threshold parameters that maximize long-term energy production. An attempt has 
been made to generalize these results from local site and turbine characteristics to 
obtain general guidelines for control algorithm design. 


INTRODUCTION 


A computer model has been developed using 
real data to simulate power production 
of a vertical axis wind turbine (VAWT). 

The purpose of the model is twofold. 

First, to develop guidelines in control 
algoritTim strategies. This basic step 
was Initially undertaken with using data 
particular to the Sandia/DOE T7-m VAWT. 

These Initial results are very Interest- 
inc and have been useful in the automa- 
tic control program for the 17-m. 

The second purpose of the model is to 
provide more general i nformati on on the 
e7 er gy a va 1 1 a bl e to a VAWT that is lost 
in the automatic control process. This 
information will, in turn, be used in 
the economic selection process of VAWTs 
for site specific cases. 

The automatic control simulation model is 
discussed In section 1, the algorithms 
on the Sandia/DOE 17-m turbine are de- 
scribed in section 2, and thefirst step 
results are presented in section 3. In 
section 4, the specific first step re- 
sults are generalized in an attempt to 
be applicable In a turbine selection pro- 
cess . 

THE COMPUTER MODEL 

The computer model was ori ginal ly de- 
signed to simulate the 17-m VAWT in auto- 
matic control in Albuquerque, New Mexico, 
and as such was based on two sets of ex- 
perimental data. The first data set is 
the electric power versus windspeed data 
collected using the method of bins CRef, 1) 
for the 17-m operating at 50.6 rpm. A 
plot of the power curve used appears in 
Fig. 1 . 

The second set of data used is a collec- 
tion of 890 hours of .5 Hz time series 
data of windspeed over 43 periods recor- 


ded during weekends and nights of Summer- 
Fall 1979. The probability density func- 
tion is plotted versus windspeed for the 
wind data in Fig. 2. The mean value of 
windspeed for the wind data is 11.1 mph 
which is above the actual Albuquerque 
mean annual windspeed. 

The actual computer model is a Fortran 
program in which the electric power versus 
windspeed data are used to determine the 
power available at every point of wind 
data. If the turbine is stopped, the 
control algorithm will use tliis informa- 
tion to decide if the turbine should be 
started. If the turbine Is running, the 
power is summed to determine the energy 
generated . 

The overall energy output that is compu- 
ted for a particular control algorithm de- 
pends on how long and when the turbine 
Is running, as well as the number of 
starts. The number of starts is used 
by subtracting the energy consumed in 
all the starts from the calculated energy 
produced when the turbine is running. For 
the 17-m turbine, each start consumes 
approximately .6 kW-hr. 

The output of the auto control -simulator 
1$ given as an algorithm efficiency, and 
a number of starts occurring over the 
duration of the wind data. The notion of 
algorithm efficiency will be given a pre- 
cise meaning following the preliminary 
notion of energy available to a wind tur- 
bine. 

Energy Available. The energy available 
to a wind turbine is that amount of energy 
that would be produced by a turbine operat- 
ing under an Ideal algorithm. This quan- 
tity may be calculated from the electric 
power data directly by Integrating the 
power produced over all times the power 
is positive, or by using the windspeed 
distribution function for positive powers. 
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Algorithm Efficiency . The algorithm effi- 
ciency is the percent the energy produced 
is of the energy available for a turbine 
operating under the algorithm control. 

ALGORITHMS TESTED 

Five different algorithms have been tested 
in parametric studies by the automatic con- 
trol simulator. The first four algorithms 
have two parameters which are adjusted to 
find the optimum operating values. These 
two parameters are referred to frequently 
and will thus^ be formally defined. 

Turn On/Off Threshold . The turn on/off 
threshold is the minimum value of average 
windspeed or average power that is required 
to start/stop the turbine by a decision 
loop of a control algorithm. 

Test Window . The test window is the dura- 
tion of the test period over which the 
average windspeed or power is calculated 
for comparison with the turn on threshold. 

We may now name and describe the five algo- 
rithms tested as follows. 

V) Discrete Windspeed Averages - At the 
end of each test window, the average wind- 
speed over the test window is compared to 
the turn on threshold. If the test fails, 
a new average Is begun and not tested again 
until the end of another test window. If 
the average windspeed is above the turn 
on threshold, the turbine is started. Dur- 
ing turbine operation, the average windspeed 
is tested against the turn off threshold 
at the end of each test window for a stop 
deci si on . 

2) Moving Windspeed Algorithm - After 
each point of windspeed data, a new average 
windspeed is calculated over the most re- 
cent test window. If the turbine Is off, 
the average windspeed is tested against the 
turn on threshold to decide on a turbine 
start. If the turbine is on, the average 
is compared to the turn off threshold. 

After a start or stop, a 60 second blind 
period is introduced to simulate an actual 
start or stop. 

3) Moving Power Algorithm - At each wind- 
speed point, a value of power is read from 
the electric power versus windspeed data. 
From the power data, a running power aver- 
age is calculated over the test window and 
compared to the turn or on turn off thres- 
hold for a start or stop. 

4) Discrete, Double Power Algorithm - At 
the end of each test length, two power 
averages are calculated. One for the 
average over the entire test window, and 
second, for the last one-tenth of the test 
window. Both averages must be greater than 
the threshold y a lue fo r a start to occur . 
Stops are based on a test of the average 
power over the entire test window. 


5) Canadian Coast Algorithm - Certain 
Canadian VAWT systems have utilized a 
semi -mechan 1 cal control system with an 
overrunning clutch on the high speed 
shaft to permit the rotor to coast when 
below synchronous speed. This system 
could conceivably reduce motoring losses 
in low wind conditions and simplifies the 
turn off condition which can be based 
on a simple rpm measurement. 

From a cold stop, one of the first four 
algorithms is used to decide if a start 
should be initiated by a starting motor 
which takes the turbine to some fraction 
of the synchronous rpm. The turbine is 
then allowed to coast until either synchro- 
nous speed is reached, at which time the 
generator engages, or 5 minutes has 
elapsed and the turbine falls below a 
turn off threshold value of rpm. If 
synchronization is reached, when the 
torque at the high speed shaft becomes 
negative, the generator disengages and 
the coast cycle begins again. Whenever 
the rotor is coasting, the system genera- 
tor is assumed to motor without shaft 
load at synchronous speed. Windage los- 
ses in the generator are accumulated 
during the coast period. The only para- 
meter is the turn off threshold rpm. 

RESULTS OF PARAMETRIC STUDIES 

The results of the parametric computer 
runs for the five algorithms are plotted 
in Figs. 3 through 7. The maximum curves 
from each of the algorithms are plotted 
together in Fig. 8 to facilitate compar- 
ing the results. The following points 
are clear from the figures. 

1) For each algorithm, a choice of the 
parameters can be made which will maxi- 
mize the algorithm efficiency, and mini- 
mize the number of starts or coast time. 

2) The algorithms may be ranked accord- 
ing to the maximum algorithm efficiency 
as follows: 

a) Canadian coast algorithm - 93.7% 

b) Moving power algorithm - 93.4% 

c) Discrete, double power algorithm - 
92.8% 

d) Discrete windspeed algorithm - 
92.3% 

e) Moving windspeed algorithm - 91.6% 

3) A bad choice of parameters, particular- 
ly the test window, will greatly reduce 
the algorithm efficiency. 

APPLICATIONS TO VAWT ECONOMICS 

The algorithm control study of the pre- 
vious sections is useful In maximizing 
energy production once a particular tur- 
bine and site have been selected. How- 
ever, if a turbine is badly matched to a 
site, the maximum algorithm efficiency 
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i ~~ ~~~ ~ ~ ~ niay be as low as 72% as shown in an exam- 
ple below. It Is clear then that control 
algorithm considerations will have an ef- P(V) 

feet on the economics of wind energy produc- 
tion. 


-6kW for 0 <V<Vq-2.2 

2.7*(V-Vq) for V^-2.2<V<3 *Vq-6,6 

5.4*{Vq- 18) for 3*V^-6.6<V. 


In order to utilize control algorithm pro- 
p^rtles in wind energy economics, it is 
nac|'l$|ry to determine if the maximum algo- 
rTthm efficiency depends functionally upon 
a suitable composition of a local windspeed 
probability distribution function (PDF) 
with a turbine electric power curve (EPC), 
Since the windspeed PDF and the turbine EPC 
are both functions of windspeed, It is 
natural to construct the energy distribu- 
tion function (EOF) from the two, which at 
any v/indspeed gives the fraction of the 
annual energy produced by the machine at 
that wjndjpeed. The EOF i $ constructed 
as tTTe normalized product of the PDF and 
the'fPC, We will concentrate then on^ 
pGah t1 fy i ng the dependence of the maxi- 
mum algorithm efficiency upon the EDF. 

The five types of algorithms discussed 
in section 2 were initially tested using 
the EPC from the Sandia/DOE 17-m VAWT. 
the EPC was determined experimentally by 
the method of bins, and the windspeed time 
series was recorded at the 17-m test site 
in Albuquerque, New Mexico. Clearly, a 
generally applicable relation between the 
maximum algorithm efficiency and the EDF 
cannot be found from the results of one 
specific set of turbine and windspeed time 
series data. Therefore, an expanded 
data base has been acquired including one 
additional set of windspeed time series 
data, and several other sets of turbine 
EPCJs, The windspeed time series consists 
of a year of data recorded by the National 
Severe Storm Laboratory from the WKY tower 
in Oklahoma City, Oklahoma from October 
1, 1976 to November 26, 1977, Most of the 
readings were taken at a sample rate of 0,1 
Hertz. The data set includes windspeed, 
w1 nd d 1 recti on , temperature, and vertical 
velocity at 7 elevations, as well as pres- 
sure,. wet bulb temperature, pyronometer , 
and rain gauge readings from selected 
elevations. The data were transcribed to 
9 track magnetic tape and made available 
for dissemination by the National Center 
for Atmospheric Research. The windspeed 
data recorded at an elevation of 26 meters 
were used In the model since this eleva- 
tion closest matches the equatorial eleva- 
tion of VAWTs of interest. The PDF of 
the windspeed time series is closely ap- 
proximated by the Rayleigh distribution 
with a mean windspeed of 10,4 mph (Fig. 

9 ). 

In addition to the Sandia/DOE 17-m EPC, 
two other types of EPC’s have been used. 
The first type is generated by the com- 
puter code PARE? (Ref. 2) for a 17-m 
V/\yx operating at 25, 30, 35, 40, and 45 
rpm's. The second set of EPC’s come from 
a one parameter family of piecewise linear 
curves defined by the function 


The parameter Vq is the cut in windspeed, 
satisfying P(Vq; = 0.0. This family of 
EPC*s was chosen to give a wide varia- 
tion In turbine characteristics depend- 
ing on only one parameter. 

Automatic control simulator tests con- 
ducted subsequent to the Introduction 
of the additional windspeed time series 
and turbine EPC's have been made using 
only the "moving power" algorithm since 
this algorithm so clearly outperformed 
the three other purely electronic con- 
trol algorithms, and was nearly equal 
to the coast type algorithm. 

The dependence of the maximum algorithm 
efficiency on the EDF may take any num- 
ber of possible forms. After finding the 
maximum algorithm efficiencies for several 
cases by running many tests of the auto- 
matic control simulator corresponding 
to different data sets, plots were made 
of the maximum algorithm efficiency 
versus various properties of the EDF’s. 

Two clear relationships emerged in the 
plots, either of which may be used to 
determine the approximate maximum algo- 
rithm efficiency corresponding to a^ 
given EDF, The first relationship is 
between the maximum algorithm efficiency 
and the slope of the EDF at the cut in 
windspeed, denoted Eq. The plot of the 
maximum algorithm efficiency versus Eq 
appears in Fig. 10. The second rela- 
tionship 1s between the maximum algorithm 
efficiency and the ratio 

V - 

R = ^ , ( 1 ) 

where Vni is the windspeed at the maximum 
EDF, and Vo> as before. Is the cut in 
windspeed. The plot of the maximum 
algorithm efficiency versus R appears 
in Fig. 11. 

From Figs. 10 and 11 it is apparent that 
the maximum algorithm efficiency may be 
given as a function of either Eq» the 
slope of the EDF at Vq> or the ratio, R, 
defined in Eg. (1). This functional 
relation may be used in a VAWT economic 
selection process for a given site, or 
more basically, in determining if a VAWT 
is badly suited controlwlse for a tur- 
bine site, as the Sandia 17-m turbine 
is for Oklahoma City with a maximum 
algorithm efficiency of 72%. 
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FIGURE 7 


CW^PARISON Of THE BEST CURVES 
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FIGURE 9 - Comparison of WKY Tower PDF with Rayleigh PDF 
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